Fragmentation with a Steady Source 
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We investigate fragmentation processes with a steady input of fragments. We find that the size 
distribution approaches a stationary form which exhibits a power law divergence in the small size 
limit, P OCI {x) ~ x~ 3 . This algebraic behavior is robust as it is independent of the details of the 
input as well as the spatial dimension. The full time dependent behavior is obtained analytically 
for arbitrary inputs, and is found to exhibit a universal scaling behavior. 
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Fragmentation underlies numerous natural phenomena 
The quantity being "split" can be the mass, mo- 
mentum, or the area, and typically, fragments continue 
splitting independently of each other. Examples include 
olymer degradation ra, energy cascades in turbulence 
, breakup of liquid droplet s B a nd atomic nuclei jIo| , 
martensitic transformations ]11| , [12| , shattering of solid 
objects HH, and meteor impacts. Fragmentation also 



arises in several topics of computer science |15j |17| . 

The simplest fragmentation models assume that the 
rate by which fragments are produced is a function of 
their size only [p~8|— p2[| . In this study, we focus on the 
classic "random scission model" where the cutting is uni- 
form and hence a fragment is cut with a rate proportional 
to its size. In particular, we are interested in situations 
where the system is subject to a steady input of frag- 
ments. Such "open" systems may be realized physically 
in impact fragmentation of solid objects where it is pos- 
sible to constantly supply the system with unfragmented 
objects. 

Fragmentation in open systems has received much less 
attention than fragmentation in closed systems. We will 
show, however, that fragmentation with input is actu- 
ally simpler than the classical counterpart as the system 
reaches a stationary state which is remarkably robust. 
Specifically, fragmentation with a steady source is char- 
acterized by an algebraic divergence of the size distribu- 
tion in the small size limit, and this behavior is indepen- 
dent of the particular form of the input. Additionally, 
the time dependent behavior, obtained analytically for 
arbitrary inputs, follows a scaling behavior. These two 
features are shown to be closely related. 

We start with a one-dimensional fragmentation pro- 
cess subject to constant input of segments. Here "one- 
dimensional" means that the fragments are characterized 
by a single variable which we shall call "length" (the frag- 
ments can be viewed as segments). Let the system be 
initially empty and intervals whose length is within the 
range (x, x + dx) are added with rate f(x) dx. Addition- 
ally, intervals are cut with a constant spatially homoge- 
neous rate; we set this rate equal to unity without loss of 
generality. Fragmentation with input has a natural geo- 
metric interpretation. Consider the segments as part of 
an infinite line. The fragmentation process is equivalent 



to deposition of point "cracks" on the line. The line is 
initially "immune" to fragmentation, but then segments 
of length x become "susceptible" to fragmentation with 
rate f(x). Hence, fragmentation with input is equivalent 
to inhomogeneous fragmentation of a growing line. 

The density P(x,t) of intervals of length x at time t 
evolves according to the following rate equation 



dP(x,t) 
dt 



-xP(x,t) 



dyP(y,t) + f(x), (1) 



The negative term on the right-hand side accounts for 
loss due to fragmentation with the rate equal to the frag- 
ment size since the cutting is uniform. The gain term 
gives the increase in fragments of size x due to cutting 
of longer fragments. The last term accounts for input of 
fragments of size x. 

The size distribution can be determined by applying 
the Mcllin transformation. The Mcllin transform (or 
moment) of the distribution, M(s,t) — J dx x s ~ 1 P(x, t), 
satisfies 



dM(s,t) 
dt 



-M(*+l,t) +/(*), 



(2) 



where f(s) = J dx x s ~ 1 f(x) is the Mellin transform of 
the input density f(x). Although this hierarchy of equa- 
tions is infinite, its linear nature makes it tractable, as 
will be seen below. 

We first examine what happens when t — > oo. In this 
limit, the length density should approach the stationary 
distribution, P(x,t) — > P 00 (x). Setting the time deriva- 
tive in Eq. (J2j) to zero gives the corresponding transform 

Moo(s) = (l + ^\ f(s - 1). Note that f(s + n) is the 

Mcllin transform of x n f(x), and (s — m)^ 1 f(s) is the 
transform of x~ m J dy y m_1 / (y) . These two facts allow 
to perform the inverse Mellin transform and yield Poc(x) 
explicitly in terms of the input function 



Poo(x) =x- 1 f(x) + 2x- 



dyyf(y)- 



(3) 



In the small size limit, the integral on the right-hand side 
of Eq. (||) approaches the average length added per unit 



1 



time, A = /(2) = J dxxf(x). Thus, the length density 
becomes purely algebraic 



Poo (2) -» 2\x~ 



when 



0. 



(4) 



This behavior is robust as the first term on the right- 
hand side of Eq. (||) always diverges slower than x~~ 3 in 
the limit x —> (otherwise, the total length input rate, 
/ dxxf(x), would be infinite). For a class of input den- 
sities, the algebraic behavior may not be necessarily re- 
stricted to small sizes. For example, for monodisperse 
inputs f(x) = XS(x — 1), the algebraic behavior extends 
to all sizes x < 1, Poo (a?) = \x~ 1 S(x — 1) + 2Aa; -3 . 

The general algebraic behavior should be contrasted 
with the exponential length distribution found generally 
in the absence of input. Algebraic distributions have been 
observed experimentally in fragmentation of solid objects 
such as rods, spheres, and bricks |T^|l4|]. Although the 
corresponding exponents measured in these experiments 
are significantly lower, typically between 1 and 2, it is 
worth noting that a steady source of fragments can serve 
as a mechanism for generating algebraic distributions. 
Curiously, algebraic distributions with an exponent close 
to 3 were reported recently in social systems which can 
be viewed as open ones (distributions of citations, of the 
number of links to sites on the internet, etc.; see e.g. 
Refs. HHH). 

The limiting size distribution is ultimately related to 
the time dependent behavior. This can be demonstrated 
using the following heuristic argument. From Eq. Q), 
the total length L(t) = M(2,t) grows linearly with time 
Lit) — A, and hence, L(t) = At. Similarly, the total num- 
ber of fragments, Nit) — M(l, f) satisfies N(t) = \t + fi, 
where /j, = /(l) = / dxf[x) is the number of segments 
added per unit time, and consequently, N{t) — ^\t 2 + fit. 
These two time dependent results imply that the average 
fragment length, (x) = L/N, decreases with time accord- 
ing to (x) ~ f _1 . For the length distribution to follow 
a scaling form, the corresponding scaling variable must 
be x/(x). The prefactor is fixed by the total number of 
fragments, N(t) ~ At 2 , so the scaling form reads 



P{x,t) ~ Xt 3 F(xt). 



(5) 



This scaling form would be consistent with a time in- 
dependent limiting distribution only when P(£) ~ £~ 3 , 
thereby implying the algebraic divergence (^). 

The full time dependent solution can be found using 
the Charlesby method This method starts with a 

formal expansion of the Mellin transform, 



t" 



M(s,t) = y2-M n (s), 



(6) 



71 = 



and proceeds by solving for the functions M n (s) itera- 
tively. Indeed, substituting the expansion (||) into Eq. (||) 
and equating similar powers of time yields Mq(s) = 0, 
Mi(s) = f(s) and M n+1 (s) = -2=2M„(« + 1) for n > 2. 
Solving this set of equations recursively gives 



M n+1 {s) - (-1)" 



(*-!)(* -2) 
(s + n- l)(s + n - 2) 



f{s + n). 



To take advantage of the inversion rules used to obtain 
Eq. (H), we re- write M n+1 (s) as 



M n+1 (s) = (-1)" 



n(n + 1) n(n — 1) 



•71-1 



f(s + n). 



From Eq. (^|), the size distribution can be written as a 
power series 



P(x,t) 



^ (n + 1) 

n=0 V ; 



(7) 



where the inverse transform of M n+ i(s) has been conve- 
niently written as (— x) n P n (x). The three terms in the 
above expression for M n+ i{s) can be inverted using the 
rules outlined before Eq. (0). The final expression for 
P n (x) reads 

Pn(*)=A(a:) + =^/ a (*) + =^M*), (8) 



with 



fi(x) = f{x), 
Mx)= dyfiy), 



(9) 



h(x) 



dyyf(y)- 



Summing the three terms separately gives the fragment 
size distribution 



POM) 



3 

J2t k fk(x)F k (xt), 

k=l 



(10) 



with the scaling functions 



F 1 {z) = z-\l 
F 2 (z) = e-', 



(11) 



Fs(z) 



[2- (2 + 2z + z 2 )e" z ] 



The function F^(z) has been obtained from the power se- 
ries F 3 (z) = J2n>0 „!(„+3) " 0ne Can Veri fy that tll<3 P re " 

vious results for N(t), Lit), and Poo(aO agree with this 
solution. Thus, we have obtained the full time dependent 
solution for an arbitrary time independent input fix). 

The size distribution of Eq. (|l(]) exhibits scaling. In- 
deed, in the limit t — > oo, x ~ > with the scaling vari- 
able z = xt kept finite, the third term in the sum on 
the right-hand side of Eq. (|o|) dominates, and the an- 
ticipated scaling behavior of Eq. (^) is confirmed with 
P(z) = F^iz). Interestingly, the only parameter relevant 
asymptotically is the overall length input rate A. 



2 



The limiting behaviors of the scaling distribution are 



F{z) 



i _ 1 

3 4 

2z~ 3 



z z < 1, 
Z > 1. 



(12) 



In particular, the large z behavior implies the correct 
asymptotic Poc(x) ~ x~ 3 , in agreement with Eq. (Q). 
Thus, for sufficiently large fragments, X 3> t , the distri- 
bution has already reached the final limiting form, while 
smaller sizes are still created. 

The formal solution ( |l0| ) has an interesting "staircase" 
structure, a time power series whose terms are products 
of time independent functions fk (x) and time dependent 
functions Fk(xt). In fact, the solution for the random 
scission model in the absence of input is also character- 
ized by a similar structure. Indeed, consider the evolution 
equation 



dP(x, t) 
di 



= -xP(x,t) + 2 dyP(y,t) 



(13) 



corresponding to the above fragmentation process in 
the absence of input. Given the initial conditions 
P(x,0) = f(x), the solution can be obtained following 
the same steps that led to Eq. (10). Again, the full time 
dependent solution is a three term expansion: 



P(x,t) 



3 
k=l 



1 fk{x) F k {xt). 



(14) 



The time independent functions are given by the same 
expressions (|]) as in the input case (with fk replaced 
by fk), while the time dependent functions are different 
from @: F x (z) = F 3 (z) = e~ 2 , and F 2 (z) = (2 - z)e~ z . 
In the limit t — > oo, x — » with z = xt kept fi- 
nite, the scaling behavior emerges again. Specifically, 
P(x, t) ~ Xt 2 F(z) with A = J dxxf(x) and the exponen- 
tial scaling function F(z) — e~ z . 

The above solutions for the input case with empty ini- 
tial conditions and no input case with arbitrary initial 
conditions can be used to construct the general solution 
for Eq. (|l|). Indeed, the sum of the solutions ( |Io| ) and 
(fli]), P(x,t) + P(x,t), is the solution for a fragmenta- 
tion process with input f{x) starting from an initial dis- 
tribution f(x). As expected, the initial conditions are 
"forgotten" in the long time limit as P(x, t) given by 
Eq. ( |To| ) dominates over P(x,t) given by Eq. ([l4|). In 
particular, the scaling solution (JsJ) is recovered, and the 
Poo (a;) ~ 2Ax -3 divergence of the limiting distribution 
holds in general. 

To examine the robustness of the algebraic behavior 
above, we consider a natural generalization to d spa- 
tial dimensions |p6|-p^|. Given that the most interest- 
ing long time behavior is independent of the details of 
the source term, we focus on the simplest monodisperse 
inputs, namely, unit hypercubes. For instance, in two di- 
mensions we add unit squares with rate A. A unit square 



is divided by choosing a point (0:1,2:2) with a uniform 
probability density, and cutting the original square into 
four rectangles of sizes x\ x x 2 , x\ x (1 — x%), (1— x%) XX2, 
and (1 — x%) x (1 — x 2 )- Similarly, the process is repeated 
with rectangular fragments. 

Let P(pc,t) with x = (x%, . . . ,Xd) be the distribution 
of fragments of size x% X • • • X xj at time t. This quantity 
evolves according to the rate equation 



d_ 

at 



|x|) P(x, t) = 2 d J dy P(y, t) + A5(x - 1). (15) 



Here, we used the shorthand notations 1 = (1,...,1) 
and |x| = x\ ■ ■ ■ Xd- The d-dimensional Mellin trans- 
form, M(s,t) = fd-xx s ^---x d d ~ 1 P(x 1 t) with s = 
(si, . . . , Sd), reduces Eq. ( |15| ) to 

dM(a,t) = { 2 d - Sl ---s d \ 

at \ s!---s d J 



M(s + !,<) + A. (16) 



We focus on the limiting size distribution P^ he) . 
Its Mellin transform M 0O (s) is found from Eq. ( Jlq ) 
by setting the time derivative to zero. One gets 

Afoo(s) = A fl + (si _ 1) ...( Sd _ 1) _ 2 d N ) • Inverting this rela- 
tion yields llf 



Poo(x) = A [<5(x - 1) + 2 d |x|- 1 $ d (0] 
with the shorthand notations 



MO 



E 

n=0 



and £ = 2 



d 

n^4 



i/d 



(17) 



(18) 



In one dimension, $i(£) = = x~ 2 , and we recover 
the one-dimensional result Poo{x) = 2Xx~ 3 . In two di- 
mensions, $2(0 = Io(2£,) where Iq is the modified Bessel 
function, and in general 3>d(£) can be expressed in terms 
of hypergeometric functions. 

The small size behavior of Poo(x) can be obtained by 
using the steepest decent method. The leading tail be- 
havior, $(z(£) — (27t£)t~ e^ d for £ S> 1, corresponds to 
the case when at least one of the lengths is small, i.e., 
ii < 1. Returning to the original variables, we re- write 
the above asymptotic as 



a/d 



P oc (x) ~ Ixj-^lnxl— aa-exp 2d(|lnx|) 



where |lnx| = Ilj=i m 7"- Thus, the fragment distribu- 
tion exhibits a "log-stretched-exponential" behavior. 

Let us consider the limiting volume distribution 
Poo (V) defined via P x (V) = J dxP^x) 6 {V - x x ■ ■ ■ x d ). 
Its Mellin transform, Moo(s) = / dV V s ^ 1 P oc (V), imme- 



diately follows from Eq. ( [16 
Using the identity (a d — 1) 



Moo(s) = A 



(5=1? 



where C 



^ x Etk fc («-c fc )- 



is the primitive d th root of unity, we 



can express Moo(s) as a sum over simple poles at l + 2£ . 



3 



Consequently, the inverse Mellin transform is given by a 
linear combination of d power laws, 
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Poo{V) = A 



d-l 



-1-2C 



fc=0 



(19) 



One can verify that the volume distribution is real since 
it equals its complex conjugate. The small-volume tail of 
the distribution can be obtained by noting that the sum 
in Eq. (JL9|) is dominated by the first term in the series, 



for V -> 0. 



(20) 



Thus, the same V -3 algebraic behavior occurs in all spa- 
tial dimensions. Clearly, this divergence is general. In- 
deed, the time dependent evolution equation ( |i~5| ) implies 
that the overall volume grows linearly, and that the over- 
all number of fragments grows quadratically. Therefore, 
the heuristic scaling argument leading to Eq. (||) extends 
to higher dimensions, and consequently, the limiting be- 
havior is given by Eq. (pC|). 

In summary, we have studied random fragmentation 
processes in the presence of a steady source. We have 
solved for the full time dependent behavior in terms of 
the input function. In the long time limit, the size dis- 
tribution exhibits a universal scaling behavior. The lim- 
iting distribution diverges algebraically according to a;" 3 
in the small size limit. This behavior is robust. It applies 
to arbitrary inputs, and it extends to higher dimensions 
as well. Interestingly, the only asymptotically relevant 
parameter is the total volume added per unit time. Ad- 
ditionally, we have shown that the scaling behavior can 
be used to predict the algebraic nature of the final size 
distribution. Hence, the scaling behavior and the limit- 
ing distribution are closely related. 

The solution for the time dependent behavior exhibits 
an interesting staircase structure. The two progressively 
weaker corrections to the leading behavior are of the 
order t~ x and t~ 2 , respectively. Such staircase struc- 
tures may be a useful tool for treating similar integro- 
differential equations which are expected to exhibit a 
scaling asymptotic behavior. For instance, one can check 
whether substituting such an ansatz leads to a closed 
system of equations for the time dependent and time in- 
dependent functions. 

The above results can be extended in a number of 
ways. One may try to derive a general solution for P(x, t) 
in higher dimensions for arbitrary input rate /(x). We 
anticipate that geometric features of the fragments will 
be interesting. Indeed, in the no input case the volume 
distribution exhibits an ordinary scaling behavior while 
multiscaling asymptotic behavior underlies the full mul- 
tivariate size distribution. Additionally, a nontrivial set 
of conservation laws exists as all moments M(s*,t) with 
n<=i s t — 2 d are conserved, as seen from Eq. (^6[). In 
the presence of input, the same moments grow linearly 
in time but multiscaling should still hold asymptotically. 
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